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Abstract 

We present Vlasov's equation and its association with Poisson's equa- 
tion in the context of modelhng self-gravitating systems such as Glob- 
ular Clusters or Galaxies. We first review the classical hypotheses of 
the model. We continue with a presentation of the Hamilton-Morrison 
structure of Vlasov's equation to study the equilibrium and the stability 
of self-gravitating systems. Finally, we present some preliminary results 
concerning some properties of the time dependent solutions of the Vlasov- 
Poisson system. 

1 Hypotheses 



It is well known that Vlasov's equation can be used to describe the dynamics of 
a self gravitating system. We recall the usual descriptive hypotheses. 

First, we consider a set of N point mass particles of individual masses 
Tni=i,..,N and we consider the associated N particle distribution function in 
the phase space /^^^ such that, at any time t 

N 

1= //W(ri,... ,r^,i) n^r. , (i) 

where Ti = (Pi,qi) is the 6- vector which indicates tlie position and the 
impulsion pi of the i*'* particle. Provided that cacli pi and are Hamiltonian 
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conjugated variables and all the forces arc conservative, one can obtain directly 
from this normalization equation a conservation equation (Liouville's theorem) 

dt ^ \m- dcii dcii- dp J ^ ' 

where Ui is the potential energy at the origin of the force fj seen by the i*^ 
particle 

dpi _ _dUi 

dt~ 9q, ^-^^ 
If the only force considered is gravitation, we have 

JV 

In order to obtain a gravitational Vlasov's equation, we first need to introduce 
marginal distribution functions /'^^and /'^•' respectively for one and two par- 
ticles as 

N 



/(I) (r,, t)= I /(^) (Fi, . . • , Tjv, t) n dT^ (5) 

N 

/(2)(ri,r,-,t)= / /W(ri,... ,r;v,t) H (6) 

We then make use of the Boltzmann ansatz which assumes the existence of a 
two-particle correlation function k (FijEj,*) such that 

Vij^j = 1, .., N {Ti, rj,t) = /(I) {Ti,t) /(I) (r,-, i) + k {Ti,Tj, t) (7) 

We now have suppose that the particles are indiscernible^, and have the same 
probability law / (F, t) = /'^^ (F, t) /N. Some direct calculations then give 



df p df d^/j df ^ f dk{T,T',t) (q-qQ 
at m aq aq ap J op |q — q T 



where ip (q, t) is the mean field gravitational potential defined from / by the 
inverse Poisson equation 

^|J{cl,t) = -Gm^l^^dT' (9) 

It is well known (see e.g. Binney and Tremaine[12]) that in the case of self- 
gravitating systems the time Tr taken by collisions to relax it, can be evaluated 
by the Chandrasekhar formula 

N N 1/9 

rT<i - ^ {Gp)-"' (10) 



Jn-ZV InAT 



^In particular this means that all particles have the same mass m. 
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where Td is the dynamical or crossing time of the system and p its mean mass 
density. The coUisional relaxing time then ranges from lO^'^^yr for Globular 
Clusters {N ~ 10^) or lO^'^^yr for a typical galaxy [N ~ 10^°) to lO^^'^'^yr 
for Galaxy Clusters (iV ^ 10^) . Hence, during a reasonable time of about a 
few hundreds dynamical times one can actually consider that two-particle cor- 
relations (namely gravitational collisions) are dynamically unimportant, and by 
consequence that 

fc(r,r',t) = o (11) 

During this long period the mean field dynamics of a self-gravitating system are 
governed by the Vlasov-Poisson system 

df p Of dib df ^ 
at m aq aq ap 



^{ci,t) = -GmJ^^^^dr' (13) 



2 Vlasov-Poisson as a Hamiltonian system 

It is often quoted that Vlasov's equation is a Hamiltonian system. It is true but 
the formulation of such a problem is generally unclear. Using classical Poisson 
brackets defined for any scalar fields A (q, p, t) and B (q, p, t) by 

(A B} ■= — —- — — (14) 
' 5q dp dq dp 

one can write Vlasov's equation in the form 

^+ P ^_^^^ = ^ + {f m = (15) 
dt m '(9q dq dp dt ' 

where 

2m 2m J |q-q'| 

is the mean field energy per particle. It is a fact that relation (15) is not a Hamil- 
tonian equation. As a matter of fact, general Hamiltonian equations give the 
total derivative of a canonical variable from a lie bracket between this variable 
and a propagator. Nevertheless, a series of works initiated by Morrison[l] , have 
shown that a Hamiltonian formulation of non dissipative physics is generally 
possible. Let us recall these results in the context of Vlasov's equation. 

We consider any physical functional of the distribution function, in a very 
general context this functional can be written 

P[n = j ^(/,r)(ir (17) 

where is a very general suitable function, for example if if do not depend on 
r, F[f] is called a Casimir and it can be easily proven that this is a conserved 



E = 
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quantity during the Vlasov dynamics. By direct calculation we have 

dF _dF _ [dfd^ 
'dt--dt-J 'dtdf'^^ ■ 

It is important to note here that 

^ = ^ (19) 

where 6/Sf indicates the usual functional derivative defined in the case of (17) 
by the relation 



SF 

—6f + 0{Sf) 



dT (20) 



F[f + Sf]-F[f]= J 
Using (18), (19) and (15) , one can then write 

Using definition (20) , one can verify that the total energy contained in the 
system 

which is a functional of the form (17), is such that 

Then, integrating the RHS of (21) by parts^ and employing (23) one gets 



dt J [Sf Sf 

which is the result introduced, perhaps more directly, by Morrison [1] . This last 
equation can be formulated as 

^ = {F, H) (26) 

^This is actually a integration by parts. One can directly show indeed that for functions 
that decay sufficiently rapidly at |q| and |p| tend to oo we have 

j dTA {B, C} = ~ J dTC {B, A} = j dT {A, B} C (24) 

In the case of equation (25) surface terms vanish as the distribution function vanishes at 
infinity. 
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where we have introduced the Morrison's bracket defined for any functional A [/] 
and B [/] of the form (17), by the relation 




Some very simple calculations show that Morrison's bracket is skew-symmetric, 
verify Jacobi's identity and can be used to form a Lie algebra. Hence, equation 
(26) naturally defines the Hamiltonian structure of Vlasov's equation. We can 
note that in this formulation there is no conjugated variable to the canonical /. 



3 Some uses ... 



3.1 Generalities 



It is well known that a regular Hamiltonian formulation of a problem provides 
room for much formal analysis. Let us see what is possible with Vlasov's equa- 
tion. In a pure formal way, one can always write 



{H,.)F:=T(F) . 



(28) 



Here, the functional operator T appearing in this last equation does not depend 
explicitly on time t, because the total energy of the system H [/] is a conserved 
quantity. Hence, one can utilize the resolvant of the Ordinary Differential Equa- 
tion (28) and we have 

F[f]=cxp[-it-to){H,.)]F[f,] (29) 

= F lU] -^{H,F + ^ {H, {H, F [fo]))f=f^ 

+ O ((t - tof) (30) 

where f,, = f {T, to). It is possible to compare this expansion with the general 
Taylor expansion obtained from equation (17) which has the form 



F[f{T,t)] = I ^{f{r,t),r)dr 



-I 



^(/(r,io),r). 



1! 



df 



2! 



dp 



+ o ((i - tof) 



dT 



Hence, we arrive at 



dip if {T, to), T) 
df 

aV(/(r,t,),r) 



-{H,F[fo]) 

{H, {H,F[fo])) 



(31) 

(32) 

(33) 
(34) 
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etc ... It is well known from quantum field theory that this expansion allows 
us see that the total energy H is the time propagator of the system. In a more 
general way, Noether's Theorem says that any conserved quantity is associated 
to a symmetry of the system, hence invoking a general canonical transformation, 
the functional Vlasov's equation can be written (see Perez and Lachieze-Rey[10] 
for details) 

fHi'.if) (35) 

where the pair {K [/] , A) represents the conserved quantity (e.g. H [/]) and the 
variable related to system's symmetry (e.g. translation), respectively. 

Below, we present two natural uses of these seemingly formal developments 



3.2 Stability analysis 

We now consider a steady state solution /q (F) of the Vlasov's equation. In a 

general way. Hamilton's theory of canonical transformations says us that any 
physical perturbation of this equilibrium is generated by some functional gen- 
erator G, such that for any functional F [/] of the distribution function we have 

f = (^,G) (36) 

where e controls the perturbation (see Perez and Lachieze-Rey [10] and reference 
therein). In particular, if the perturbed state is described by a distribution 
function / (F, t) written as 

/ (F, t) = fo (F) + e/i (F, t) + e'h (F, t) + O {e') (37) 

Taking F = Id in equation (36), one can then verify that 

MT,t) = {g,fo}, f2iT,t) = ^{g,{g,fo}}, etc... (38) 

where the so-called generating function g is the functional derivative of G. This 
kind of theory says that the linearized perturbation /i := {g,fo} is generated 
from the unperturbed distribution fo (F) via a canonical transformation by some 
small perturbation g (T,t). And so on for high order non linear terms. 

Considering such kinds of perturbation is not restrictive. As a matter of 
fact, only perturbations of this form can be generated via a deformation that 
enforces all the constraints associated with conservation of phase (Liouville's 
theorem). 

Taking F = H in equation (36) allows us to obtain an energy perturbation 
induced by equilibrium perturbation 

H If] = Ho + eHi +e^H2 + (e^) (39) 
= H [fo] - e (G, H)^^^^ + (G, (G, H [/o]))/=/„ + O {e^) 

The first two terms of this expansion are very useful in stability analysis. 



6 



3.2.1 Equilibrium 

It is a well known property that if a system is in an equilibrium state, the first 
order variation of its total energy vanishes. This property can be directly shown 
from the use of the e— order term in the expansion (39). As a matter of fact, 

H, = - iG,H),^,^ = -j\^^,'^^f^dT=jf^ [e, 1} dr (40) 

an integration by parts and the use of Vlasov's equation then gives 

which is evidently a vanishing term if /„ is an equilibrium state. 
General equilibrium states are the following : 

• fo = fo (E) distribution functions represent all q— spherical and p— isotropic 

equilibria of coUisionless self-gravitating systems ; 

• fo = fo {E, LP') distribution functions where is the mean field squared 
angular momentum per particle, represent all q— spherical and p— anisotropic 
equilibria of coUisionless self-gravitating systems ; 

• fo = fo {E, Lz) distribution functions where is the mean field 2;— component 
of the angular momentum per particle, represent an equilibrium self- 
gravitating system with a 2;— symmetry in q— space. 

For detailed proofs of these affirmations and properties of these equilibria, 
we refer to Perez and Aly[8] . 

3.2.2 Equilibrium stability 

The stability of a given equilibrium described by the time-independent distri- 
bution function is classically related to the sign of the order term in the 
expansion (39) 

H2^\{G,{G,H[fo\))f^f^ (42) 

More precisely, one can prove that an equilibrium state described by a distri- 
bution function fo is guaranteed to be linearly stable if the second variation H2 
is strictly positive for all perturbations. The positivity of H2 thus provides a suf- 
ficient criterion for linear stability. If is of indeterminate sign, the situation 
is less obvious. In fact the existence of negative energy perturbations does not 
necessarily imply an instability (there exist some counterexample). However, 
one does at least expect that it guarantees non linear instability in the presence 
of dissipation. This idea goes back at least to Moser [5] . The implications for 
the Vlasov's equation were first considered by P..J. Morrison[6] in the context 
of plasma physics. That negative energy perturbations imply an instability to- 
wards dissipation was proved more recently in a fairly general context by Bloch 
et al [7] . 
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Let us indicates below the essential calculations which yields an explicit 
formula for H2. From (27) , taking g SG/df, we get 



dT 



(43) 



f = fo 

from the differential properties of the Poisson bracket we have 



Hi'') 



1 



fo{g 



<5/' 



E'\fo+\9' 



SE' 



f:, + {g',E'}Sa{T-T') 



dr'l dT 
(44) 



An integration by parts of the two first terms in the middle sum and the inte- 
gration of the dd Dirac mass term gives 



H^'^ = \j fo [g, [{g, E} + J^-^ {E', a dT' - J {g', a ^rfr' 



dT . 



If fo is an equilibrium then {E, vanishes, we remark also that from (16) 



SE' 5ip' 



|q-q'l 



(45) 
(46) 



to arrive at 



H^'' = \jfo {g, {g, E}) dr + ^Ju !^g, J ^f^dr'j dr . (47) 

A last integration by parts on each term allows us to see that 

H^'^ = -\j{9^ fo] {9, E}dT-^jj ^llM^^dT'dV (48) 

Studying of the sign of this quantity, one can study the stability of equilib- 
rium associated to fo against perturbations generated by g. Classical results 
are presented in Perez and Aly[8], Perez et al.[9] and Kandrup[2], [3], [4] . They 
can be summarized as follow : 

• If lo = fo {E) and foE ■■= dfo/dE < 0, then for all g we have F^^) > 
0. Hence, any spherical and isotropic stellar system with an outwardly 
decreasing distribution function is linearly stable against all perturbations. 

• If /o = fo {E, L') , foE < 0, and g is such that {g,L'} =0 then > g, 
such perturbations are called preserving perturbations. All radial pertur- 
bations are preserving^. Hence, any spherical and anisotropic stellar sys- 
tem with an outwardly decreasing distribution function is linearly stable 
at least against all radial perturbations. 



This proposition is not reciprocal. 
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• If /o = fo {E,L^) (rcsp. fo {E,L^)), foE < 0, and /o^s < (resp. foL, < 



0) then there always exists non radial (resp. nonaxisymetric) perturba- 
tions, corresponding to nonzero mass density perturbations, for which 
< 0. It was conjccturcd[9] that when (E, L"^) /„ (E) Sj, (L^) 
(resp. fo {E, Lz) fo (E) Sd (i^)) then if(^) < for all perturbation 
generator g. This may be associated to Radial Orbit Instability. 

3.3 Study of the transient regime of gravitational Vlcisov- 

Poisson dynamics 

The transition to a steady state is a crucial problem in self-gravitating systems 
dynamics. As a matter of fact, observations suggest that almost all galaxies 
or globular clusters manifest overall similarities and that they do not deviate 
very much from a statistical equilibrium. We may also note that numerical 
simulation of collections of point masses tend to typically evidence a rapid, 
systematic evolution, on dynamical time scale T^, towards a state, which is at 
least in a statistical sense, nearly time-independent (see Roy and Perez[13] ). 

In the toy model of one dimensional gravity, where one considers the in- 
teractions of infinite one dimensional-plane sheets, careful simulations can be 
performed (e.g. Mineau et al.[14]). The net result of such integrations is that 
typically such systems do not exhibit a pointwise evolution toward a steady state 
configuration. It is significant that solutions of ID Vlasov-Poisson system for 
the same initial conditions also exhibit an evolution toward a final state which 
is not strictly time-independent. 

In the unphysical 2D gravity model, a complete thermodynamical approach 
can be performed and show that there can exist an unique and well defined equi- 
librium state associated to a maximum entropy principle (see Aly and Perez [11] 
). This result should be soon confronted with numerical simulations of 2D 
Vlasov-Poisson system. 

In the real 3D world, one can ask if there exists a correspondence between 
time-dependent Vlasov-Poisson solutions, numerical simulations of large collec- 
tions of point masses and corresponding thermodynamical systems. 

In order to progress toward this end, let us present now, some uses of Hamil- 
tonian techniques presented above which perhaps could allow one to exhibit 
time-dependent solutions of 3D Vlasov-Poisson system. 

It is a fact that, at a given position q and for a given time t the gravitational 
potential is a functional of the distribution function 



one can then always take F — ijj in the functional Vlasov's equation to have 



hence, we can have a formal time development of the potential from which. 




(49) 



{H, .) ^ 



(50) 
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applying Laplacian operator one can deduce the temporal evolution of the mass 
density. There is no problem in principle, rather only technical ones ... 

As — {H, .) is a time-independent operator, we can employ the resolvant of 
equation (50) to write 

[/] = <l + + ^^^1 + o({t- t^f) 

+ {H, {H, V'q.t^ [fo]))f^f^ + 0[{t- tof) (51) 

The three first terms of this expansion are easily calculable. 
The first term is obtained by definition 

V-q,*,. [fo] := -Gm J ^^^^dV (52) 



From this relation one can see that 



Sf |q-q'| 
and one can compute the second term 



(53) 



= J [E\Gm\q~c^\-'}f{T',to)dT' 

= -gJ p'.^h_^f(r',to)dT' (54) 

A first hypothesis on the initial distribution function can be now posed to sim- 
plify the calculus. 

Hypothesis 1 : The initial state is symmetric in the velocity space. The 
consequence of this assumption is that for any vectorial space field B (q) and 
any vectorial odd velocity field A (p) , we have 

J A ip). Bid) f{r, to) dr = o (55) 

Under this non-restrictive assumption, one can easily show that 

<1 = (56) 
Moreover, in a general way one see also that 
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which allows the computation of the third term in the potential expansion. 

41 = (^'41),^,^ = -gJ |i^',P'.^^^^} /(r',t.)dr' (58) 

= A + B 

where 



A = -G^l4r'f(r',..)'-^.?h_^ (59, 



and 



V'q^t^ could be greatly simplified assuming some physical properties of the initial 
state / (r, to). First, we suppose that the initial density is constant on a bounded 
domain f2 C ffi.^. Then we may pose 
Hypothesis 2 



/ 



f{T,to)dp=i'o-Xn = —-Xn (61) 
m 



where Po = mvo is a strictly positive constant {uo is the value of the uniform 
initial number density whereas po represent the mass density) and xa denotes 
the characteristic function'' on f2. 

This hypothesis which can appear very restrictive, is fairly adapted for the 
description of the initial state of the physical material before it collapse to form 
a self-gravitating structure. Under this hypothesis, we have 

an integration by parts in which all surface terms vanish then gives 

A = Gpo [ ^q.,t„.A'(|q-q'r') . (63) 

It is well known that |q — q'!""*^ is proportional to the Fundamental solution of 
the 3D radial Laplacian, more precisely we have 

A'(lq-q'r') =-47r5d(q-q') (64) 

hence 

^=-%^ (65) 



^Classically, the characteristic function is defined to be unity on Q, and is zero elsewhere. 
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where we have introduced the Jeans time Tq :— {A'KGpo)~ ■ 

The second term is a httle bit more comphcated, in Cartesian coordinates 
we have 

f'it/T'/Cr'.UKpi^^^pl^ (66) 

We now make our last assumption. 

Hypothesis 3 : The initial state has an isotropic and constant velocity dis- 
persion cr^ 

Vi = l,2,3 I dp f{T,to)p'i= iyo'^^Xa = Po^Xa (67) 

Physically this assumption means that the initial state is isothermal. 
Coupling all the hypotheses, we have 

B = a^Gpo[ A'flq-q'r') =-4 (68) 

The beginning of the potential expansion may be written as 

V'q,* [/] = ^q,t„ [/o] - K,*., [/o] + CT^] + O (it - tof) (69) 

Attempting to go further in this expansion is extremely perilous. As a matter 
of fact, we see from the term A that, the functional derivatives of V'q,f„ do 
not come directly due to the /-non linearity incoming in the term i/). /. The 
calculus then becomes extremely tedious ... 

At this stage, only some general remarks can be made using such a technique 
in a general way : 

• Under hypothesis 2, all odd terms vanish in the expansion 

<!t^'^=o, VfceN (70) 

• Keeping only linear terms, it seems^ that only terms of the same nature 
of V'q^tc remains in the expansion, i.e. 

^^(2.) ^ VWo|M + , Vfc e N (71) 
neglecting the constant, this gives 

V'q,* = V'q,*. \Q E(-l)'lfc)^^ = [/o] (72) 

this maybe corroborates the r,,— periodic damped regime which appear in 
generic gravitational collapse numerical experiments. 

^It is more conjectured using preliminaj: results from high order terms, than strictly proved 



12 



Wc stress that there are prchminary results. We hope to give a more complete 
description of the transient regime of Vlasov-Poisson dynamics in an upcoming 
article. 
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